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Part  II 

Calculations  of  Steady  Diabatic  Flow 
in  One  and  Two  Dimensions 


ABSTRACT 


The  purpose  of  this  report  is  to  describe  specific  calculations  of 
diabatic  flows  in  order  that  both  the  physical  understanding  of  these 
flows  and  methods  for  their  computation  may  bo  improved. 

After  a  review  of  the  previously  formulated  theory,  one-dimens¬ 
ional  radial  and  vortex  flows  are  considered.  They  illustrate  several 
differences  betweon  adiabatic  and  diabatic  compressible  flows  such--as 
absence  of  limit  circles  and  non-minimal  stream-tubo  area  at  3onic  velo¬ 
city.  Two-dimensional  (uniplanar)  flows  describable  by  a  potential  are 
computed  next.  For  ono  form  of  potential  flow,  the  exact  partial 
differential  equation  for  the  potential  is  linear  with  constant  coeffi¬ 
cients  and  always  of  elliptic  type,  if  certain  indirect  analytical  re¬ 
strictions  are  placod  upon  the  flow  and  upon  the  variation  of  the  rate 
of  heat  addition  from  point  to  point*  Linearization  of  a  wide  variety 
of  uniplanar  potential  flows,  including  those  of  elliptic,  parabolic 
and  hyperbolic  type,  is  possible  if  tho  boat  addition  rate  is  not  too 
largo.  Scaling  laws  may  be  based  on  this,  linearization. 

A  third  type  of  diabatic  flow,  somov/hat  loss  general  than-  those 
preceding,  is  analyzod  in  more  detail.  Hero  tho  flow  is  uniform  at 
infinity  but  is  porturbed  by  a  localized  heat  sourco.  Tho  velocity, 
pressure  and  donsity  variations  aro  calculated  ovor  the  field  of  flow. 
Sinco  tho  largest  perturbation  is  -on  the  density,  it  is  the  density  per¬ 
turbation  which  limits  the  sizo  of  the  heat  source  that  still  permits  a 
first-order  perturbation  calculation. 

Th6  concluding  section  describes  a  formulation  of  the  general 
equations  for  uniplanar  flow  in  terms  of  the  components,  n,  v  of  tho 
Crocco  vector  W.  Tho  partial  differential  equation  obtainod,  say  for  u, 
is  always  hyperbolic  and  quasi-linear  if  it  is  assumed  that  v(x,y)  is 
specified  first.  Tho  heat  addition  function  is  directly  calculable  once 
u  and  v  are  known. 
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LIST  OF  SYMBOLS 


a  Telocity  of  sound 

0^  specific  heat  at  constant  pressure 

Dsl-1 2  S  +  g/ 

D  k  constant  value  of  D 

o 

e  radial  \*iit  rector 
~r 

unit  vector  orthogonal  to  9 

s(»  j  *  T/P  8  t]  1/2 


g/  =  d  In  g/d  In  N 


1,  j  unit  vectors  in  the  x  and  y  directions  respectively 

fT 

If  Mach  number 

K  vector  function  chosen  to  represent  an  irrotational  field  of 
~  flow 

n  number  of  dimensions  less  one 
n  unit  vector  normal  to  streamline 
p  pressure 

Q  heat  added  to  the  fluid  per  unit  mass  and  time 

%  (6  8  T)l/2] 

%  2  «/K  <»  -  ■*>] 

heat  added  j>pr  unit  mass,  per  velocity  squared,  per  unit  length 
>R  gas  constant 
r  radial  distance 
S  speoifio  entropy 
s^  unit  vector  tangent  to  streamline 
T  statio  temperature 

stagnation  temperature 
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stagnation  pressure 

u* ,  v'  pertubation  velocities  in  i  and  j  directions  (Section  4) 

u#  v  components  of  the  voctor  17  in  the  x  and  y  directions  respectively 
(Section  5) 


V  fluid  velocity^ 

l/2 

local  value  of  the  limiting  velocity  =  (2  c^  T^)  ' 

5^jAt 

x,  y  coordinates 

^=(/+  l)/(  7  -  1) 

* Y  ratio^of  specific  heats 

9  angular  coordinate  in  a  polar  coordinate  system 
£=!  a  function  of  u  and  v  containing  no  second  partial  derivatives 


P  density 

$  potential  for  an  irrotational  N  field 

XI  « 

*  ^ 

0^  potential  for  irrotational  V  /V^  field 

w  ■  m 

ZZ.  If  —  V 

..  vorticity  7  VMP  \/AVs  and  V*W 

1  as  a  subscripts  refer  3  to  a  transformed  coordinate  system 
(*)  refers  to  a  porturbation  of  a  variable 


INTRODUCTION 


The  study  of  diabatic  (that  is  non-adiabatic)  steady  flows  is  re¬ 
quired  in  order  that  phenomena  associated  with  combustion  in  steadily 
moving  g&ses  may  be  thoroughly  understood*  Since  the  heat  release 
during  combustion  has  a  strong  influence  upon  the  gas  flow  we  consider 
a  theoretical  model  for  moving,  burning  gases  in  which  only  the  dy¬ 
namical  effects  of  the  h  at  release  are  considered,  and  the  effects  of 
viscosity,  diffusion  and  of  change  in  specific  heat  or  in  composition 
of  the  gas  are  neglected*  The  steady  mean  flow  of  a  turbulent  burning 
gas  could,  for  axarnple,  be  described  by  this  model,  or  the  pressure 
distribution  along  a  boundary  layer  owing  to  combustion ’outside  the 
boundary  layer  could  be  computed.  Theory  based  on  the  model  is  thus 
concerned  with  steady  diabatic  flow  of  an  idealized  fluid  just  as 
classical  aerodynamic  theory  was  concerned  with  the  flow  of  air  out¬ 
side  boundary  layers.  ’ 

The  direct  problem  in  diabatic  flow  theory  consists  in  determining 
characteristics  of  the  flow  pattern  from  knowledge  of  the  heat  source 
distribution.  The  inverse  problem,  in  which  calculation  of  heat  sources 
follows  from  knowledge  of  the  flow  pattern,  is  easier  mathematically  and 
yet  is  of  sane  help  in  planning  the  solution  of  the  direct  problem.  In 
our  earlier  investigations  of  both  direct  and  inverse  problems  of  dia¬ 
batic  flow, (!»2)  fornal  manipulations  of  the  partial  differential  equa¬ 
tions  have  been  emphasized.  The  nature  of  the  resulting  physical  and 
mathematical  problems  has  been  described,  but  no  detailed  solutions  of 
the  problems  have  been  given.  Accordingly,  in  this  report  we  wish  to 
discuss  several  explicit  solutions  of  the  equations  for  steady  diabatic 
flow  and  shall  derive  from  them  a  more  intimate  understanding  of  the 
relation  between  the  physical  problems  of  aerodynamic  combustion  and 
their  mathematical  solution.  We  note  that  in  some  cases  only  part  of  a 
calculated  flow  pattern  has  direct,  physical  interest,  but  that  part  is 
worth  calculating  even  by  what  appears  to  be  artificial  methods. 

The  differences,  f$r  each  variety  cf  fluid  motion,  between  adia¬ 
batic  and  diabatic  flow  are  noteworthy.  In  adiabatic  flow  the  enthalpy 
of  each  fluid  particle  is  constant  and  the  entropy  change  of  oach 
particle  is  zero  except  across  shocks.  On  account  of  this  special 
thermodynamic  behavior,  adiabatic  compressible  flows  are  characterized 
physically,  for  example,  by  minimum  area  of  stream  tubes  at  sonic 
velocity  and  mathematically  by  the  identification  of  subsonic  and  super 
sonic  flow,  respectively,  with  elliptic  and  hyperbolic  equations.  Also 
in  this  adiabatic  case  the  two  kind?  of  compressibility  effects,  namely 
those  associated  with  high  Mach  number  and  those  associatod  with  density 
or  temperature  changes,  are  closely  coupled.  In  diabauic  flow^on  the 
tne  other  hand^  the  enthalpy  and  entropy  of  fluid  particles  can  vary 
wing  either  to  addition  of  heat  locally  by  combustion  or  by  conduction 

TTr*^;  Hicks,  "Diabetic  Plow  of  a  Compressible  Fluid”,  Quart.  App, 

Math.  6„  221-23?  (Oct.  1948)  (Referred  to  later  as  D~l). 

(2)  Lb  Hicks,  ”0n  the  Characterization  of  Fields  of  Diabatic  Flow”, 

ERL  Report  No.  633,  (May  10,  1947).  Part-  I  -  Gcnoral  Theory  of 

Steady^ Diabatic  Flow:  or  Quart.  App.  Math.  6,  407-416  (Jan  1949). 

(Referred  to  later  as  D**2). 
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of  heat  from  the  neighborhood  of  the  fluid  particles*  One  lias  there¬ 
fore  no  a  priori  idea  of  the  behavior  of  diulntti  c  flows  -with  regard  to 
minimum  area  of  stream  tubes  and  cannot  associate  the  magnitude  of  the 
Mach  number  "with  the  mathematical  type  of  the  flow.  Also  the  two  mani¬ 
festations  of  ’’compressibility”  are  now  less  tightly  coupled  since  the 
density  can  change  appr^iably  with  but  small  change  in  velocity  and  at 
small  Mach  numbers*  Rrora  the  one-dimensional  theory,  however,  we  know 
that  the  effects  of  heating  and  of  stream  tube  area  variation  are 
similar*  For  example,  a  transition  from  sub-  to  super-sonic  flow  may  be 
accomplished  in  adiabatic  flow  by  a  converging-diverging  nozzle.  In 
diabatic  flow  it  may  be  effected,  in  a  duct  of  uniform  area,  by  proper 
addition  and  abstraction  '?€  heat.  Although  one  can,  therefore,  guess  at 
some  of  the  effects  produced  by  heating  in  two-  or  three-dimensional 
flow,  it  is  possible-  to  gain  a  more  general  insight  into  diabatic  flow 
in  the  higher  number  of  dimensions  by  studying  examples  of  two-  or 
three-dimensional  theory. 

9 

In  Section  1  the  basic  equations  developed  in  D-l,  2  for  two-  and 
three-dimensional  diabatic  flow  are  summarized.  The  simplest  flows  are 
then  considered  first  (Section  2).  These  are  flows  which  are  one-di¬ 
mensional  in  the  sense  that  for  them  the  local  Mach  number  depends  on 
but  one  space  variable.  The  connection  with  both  the  one-dimonsional 
or  hydraulic  approximation  for  flow  in  ducts  and  two-  or  three-dimension¬ 
al  diabatic  flows  can  bo  exhibited  as  well  as  the  relationship  to  various 
adiabatic  flows.  The  irrotational  uniplanar  diabatic  flows  (Section  3) 
are  next  in  order  of  simplicity  because  they  are  describable  in  terms  of 
a  single  potential  function.  Special  elliptic,  parabolic  and  hyperbolic 
flows  are  discussed,  and  then  a  general  treatment  of  irrotational  flows 
is  given,  baaed  on  assumption  of  a  slightly  perturbed  flow,  which  leads 
to  similarity  laws. 

The  third  type  of  flow  discussed  (Section  4)  is  of  the  greatest 
basic  importance.  Here  one-  asks,  what  are  tho  effects  of  heat  added 
locally  in  an  unbounded  gas,  flowing  steadily,  whose  velocity  is  uni¬ 
form  far  upstream  ojp  tho  local  source  of  heat.  The  equations  are  first 
linearized  according  to  a  perturbation  scheme  since  even  for  low  Mach 
numbers  the  original  equations  are  still  non-linear.  It  is  then  possible 
to  understand  how  the  effects  of  heat  sources  distributed  in  an  arbitrary 
fashion  throughout  the  field  of  flow  could  be  built  up  by  superposition 
of  the  elementary  solutions «  The  construction  of  the  elomen'eary  solution 
itself  offers  somo  difficulty  if  tho  appearance  in  the  field  of  vortex 
filaments  and  of  infinite  changes  in  enthalpy  is  to  be  avoided. 


From  tho  consideration  of  these  detailed  problems  one  can  formulate 
what  appears  to  be  a.  reasonable  approach  to  general  two-  or  throe-  dimen¬ 
sional  problems  in  diabatic  flow.  This  formulation  has  been  examinod 
particularly  (Section  5)  for  rotational  two-dimensional  diabatic  flows,  ^ 
whore  it  amounts  to  specifying  throughout  the  field  of  flow  one  component 
of  tho  Crocco  vector  (Y l  =  V/V^  where  V^.  is  the  local  value  of  the  limit- 

velocity),  and  specifying  the  second  component  and  its  normal  derivative 
along  some  curve  in  tho  field.  The  differential  equation  to  bo  solved 
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is  then  always  hyperbolic.  After  calculation  of  the  second  component 
of  W  the  nature  of  the  heat  sources  and  of  other  characteristics  of  the 
flow  could  be  computed.  This  kind  of  procedure  in  which  the  heat  sources 
are  not  specified  first  seems  to  introduce  somewhat  more  tractable  equa¬ 
tions  ,  without  unduly  restricting  the  flow  pattern,  in  almost  all  cases 
of  diabatic  flow  th*t  have  been  studied  whether  in  one  or  more  dimensions. 
(See  Section  5  and  References  12  and  13'  for  discussion  of  new  work  along 
these  lines.) 

The  theory  developed  in  this  paper  has  been  presented  in  part  at 
meetings  of  the  American  Physical  Society  in  1947-48  and  in  the  Third 
Symposium  on  Combustion.  Flame  and  Explosion  Phenomena,  (Williams  and 
Wilkins,  1949)  (212  -  222).  The  present  report  differs  from  the  Sym¬ 
posium  account  only  in  additions  and  corrections  that  have  been  made, 
particularly  in  Sections  2  to  5. 


1.  SUiaiARY  OF  BASIC  THEORY 

For  an  inviscid  compressible  fluid  containing  heat  sources,  the 
equations  of  steady  flow  are 

y  p  +  />  V  •  VZ  =  0  (M) 

V*  f  V  =  0  (1.2) 

Op  V  -  \7Tt=TV«  VS=Q  (1.3) 

where  the  symbols  p ,  />  ,  T,  T^.  V,  S  and  Q  stand  respectively  for  pres¬ 
sure,  density,  temperature,  stagnation  temperature,  fluid  velocity, 
specific  entropy  and  hedt  added  to  the  fluid  per  unit  mass  and  time*  The 
quantity  (Q ft)  then  gives  the  boat  added  to  unit  mass  in  unit  distance 
along  a  streamline.  Tho  stagnation  temporature  is  a  measure  of  the  total 
energy  of  a  fluid  particle.  For  the  perfect  gas  here  considered  the 
speoific  heat  at  constant  pressure,  c  ,  is  constant  and  the  equation  of 
state  is  ^ 

p  -  R  f  1=0  (1.4) 

Although  variation  of  specific  heat  and  of  gas  constant  R  occur  in  com¬ 
bustion  tones  and  the  phenomena  of  diffusion  and  turbulence  often  play 
a  role,  those  complications  are  neglected  here,  as  in  earlier  diabatic 
flow  theory,  in  order  that  the  important  effects  of  heat  generation 
alone  can  be  examinod.  It  is  noted  that  the  term  Q  could  include 
effects  of  heat  conduction  explicitly  although  this  possibility  is  not 
examined  further  here. 

As  has  boen  shown  previously^0*'*  ^  the,  equations  are  placed 

in  their  most  generally  useful  form  if  transformation  is  made  to  Crocoo 
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•ectdr*  W  and  stagnation  pressure,  p^,  in  place  of  velocity  vector  and 
(static)  pressure  p  through  the  equations 

T=VtW  (1.5) 

P  =Pt,(l  -X2)^'1  (1.6) 

'■here  is  the  local  value  of  the  limiting  velocity 

Tt  =  (2  ‘p  Tt)lA  (1-7) 


and  the  stagnation  temperature  is  related  to  the  (static}  temperature 


T  =  Tt(l  -  l2) 


(1.8) 


The  transformed  equations  are 

Vlog  pt  =  [(1  -T*2)"1  Wx  (yx»)  -  W  J  (1.9) 

V  -  (1  -  W2)1^  -1  W  =  qw  (1  -  W2)1^  “1  (1  +  Zl*  w2)  (1#10) 

jf  -  SJ  log  Vt  =  (1  -  If2)  q^  (1.11) 

in  which 

%  =  «At3  (1  -  H2)  (1-18) 

*  The  Crooco  vector  was  first  used  by  the  senior  author  and  his  co¬ 
workers  at  the  Cleveland  Laboratory  of  the  MCA  in  1945.  The  vector 
|T  suggested  itself  as  a  natural  generalization  of  variables  used  in 
one-dimonsiona]  flow  theory 0f  Crocco’s  language,  and  as  a  corol¬ 
lary  of  the  Mach  vector  we  nad  introduced  in  1943'*'.  Recently  Munk 
and  others'^* ®)  have  also  introduced  the  vector  Jf  (their  "reduced 
velocity  vector")  which  they  aro  now  using  in  detailed  studies  of 
general  adiabatic  flow. 

(3) 

v  *  B.  L<>  Hicks,  D.  J.  Montgomery  and  R.  H.  Wasserman,  "On  the  One-Dimen¬ 
sional  Theory  of  Steady  Compressible  Fluid  Flow  in  Ducts  with  Friction 
and  Heat  Addition",  J.  App.  Phys.  18,  891-902  (Oct.  1947). 

^  B.  Lc  Hicks,  P.  S;  Guenther  and  R.  H.  VTasserman,  "How  Formulation  of 
the  Equations  for  Compressible  Flow",  Quart.  App.  Math.  5,  357-361, 
(Oct.  1947) 

(  5) 

1  M.  II.  Munk  and'R.  Prim*  "The  Multiplicity  of  Steady  Gas  Flows  Having 
the  Samo  Streamline  Pattern",  Proc.  Nat.  Acad.  Sci.  33,  137-141 
(May  1947). 

(gV 

M..M.  .Munk  and  R.  C.  Frim,  "On  the  Canonical  Form  of  the  Equations  of 
Steady  Motion  of  a  Perfect  Gas",  J„  App.  Phs.  19,  957-958,  (Oct.  1940). 


■t 
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If  a  flow  described  by  equations  (1.9  -  l.ll)  \7ere  irrotational  in 
the  IT  field  (therefore  admi tting  a  potential  <p  _)  then  heat  addition 
must^bo  such  that  (cf.  D-l)  * 


n,  =  q-C^)  (1.1Z) 

so  that  q^,  is  coiffetant  on  surfaces  nqrnftl  to  streamlines.  Although  this 

is  a  physically  important  and  visualizable  case,  many  other  diafcatic 
flows  rotational  in  the  77  field  nay  also  bo  derived  (cf.  p-2)  from  & 
potential  function  introduced  with  the  more  general  transformation 

V  =  ffiOO  Rt]  1/2  V#  ,  j  K  j  =  K  =  |  |  (1.14) 

in  which  g (JJ )  can  be  arbitrarily  specified.  The  accompanying  restriction 
on  Q  was  found  to  be  expressible  as 


U  -  H2C  +  if).  9,,  =jB  (Cf  -  R2£)  V  •  J 

+  g'1  a  -  ~  »2g)(i  +  \  Lf)  y  (<*u),  (i.is) 

-  d  Jcg  g/d  log  R 

qu  -  Q/4  cp  (g  R  T)1/2  (1.16) 

We  noto  that  q^  and  a„  have  always  the  same  sign  as  Q  although  their 

variation  in  the  field  of  flow  may  be  different  owing  to  the  variation 

of  the  factor  T^2  g  (K ) .  The  quantity  ^  •  jS  is  the  fractional  rate  of 
change  of  stream- tub^  area  along  a  streamline  (4),  S,  being  the  unit  vec¬ 
tor  JV/TT  tangent  to  the  streamline  at  each  point  in  the  field  of  flow. 

The  function  ?(<£,.)  controls  the  variation  of  (g  T)  along  streamlines. 


(1.17) 


J  G  j*  *  V  leg  (g  T)  -  F(0U) 

fixes  the  variation  of  pressure  throughout  the  field  of  flew, 

log  P  *  /  g  d  ~  K2  a  ^  =  constant  (l»18) 

and  enters  the  partial  differential  equation  which  tho  potential  ^ 
must  satisfy, 

U*l  L  *  '  Ur’>  i\  '  ’  U* i  iXj  MiUj 

i  H 


=  Na(N 


.2  1 

'  i  Q 

7  d*' 


^  fa*) 
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(1.19) 


The  type  or  equation  (1.19)  is  elliptic,  parabolic  or  hyperbolic  ac¬ 
cording  to  whe thcr  D  ii  =,  gi  0.  The  function  D  depends  only  on 
K  and  its  fora  is  fixed  by  choice  of  g(w)  alone. 


D=l-Sj+  £f  (1-20) 

Because  of  the  presence ^of  two  arbitrary  functions,  P(^,,)  and  g 

irrotational  t  fields  are  expected  to  be  sanewh&t  less  restricted  than 

they  appear  . to  be  at  first  sigit.  A  family  of  irrotational  fleas  para- 

•aterized  by  D  =  D  =  constant  were  described  in  D  -  2  where  exoressiens 
*  o 

sera  given  for  g(S)  for  all  rcai  values  of  5  .  It  is  converiinet  to  have 

°  12 

available  the  corresponding  expressions  for  =  exp  J  g  d  —  I<  ,  nft=*ly. 


\  =  / 1  -  (jA)1 


+  Dcy  -i 


K  t  -1 


=  -  flog  oa)]  _1 


1 


(1.21) 


Various  general  properties  of  these  equations  and  of  the  correspon¬ 
ding  physical  3ystens  have  been  discussed  in  D-l,  2.  We  shall  apply 
then  here  to  the  discussion  of  a  number  of  special  cases.  In  all  of 
these  cases,  equation  (l.ll)  shows  that,  in  the  notion  of  any  fluid 
particle,  the  iVactiodhl  rate  of  change  of  its  total  energy  with  re¬ 
spect  to  distance  is  simply  2  o  (l  -  If2 )/W  while  the  fractional  rate 

cf  change  with  respect' to  tine  is  2  (1  -  TT)  V^_. 


2.  CTJE-DI!^:SIC:iAL  FLOTS 

As  in  other  branches  of  fluid  dynamics,  an  understanding  of  one- 
dinenslonal  flows  is  essential  before  the  more  difficult  problems  of 
calculating  multi-dimensional  flows  are  attempted.  Uniform,  parallel, 
one-dimensionsional  flow  tliat  is  diabatic  has  been  the  subject  of  cany 
papers  (cf.  for  example,  reference  3  and  papers  cited  there)  because  of 
Its  direct  applicability  to  jet  propulsion  and  to  detonation  phenomena. 
We  shall  give  here  an  Introductory  discussion  of  one-dimensional  radial 
and  vortex  flows. 

Radial  Flows 
The  equation 

W  =  W(r)  c 

r 

doscribcs  a  simple  radial  flow  for  which  the  magnitude  of  the  Crocco 
vector  i7  depends  only  on  the  radial  distance  and  the  direction  of  the 
vector  Ticld  is  given  by  the  radial  unit  .vector,  c If  (n  +  l)  indi¬ 
cates  the  number  of  dimensions  (n  =1,2  correspond  to  line  and  point 
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sources  respectively)  then  the  equations  of  notion,  continuity  end  energy, 
equations  (1,9  -  1,11)  in  W  language  beeooe  simply 


^  log 


pt  -2  y 


(2.2) 


“a  rn  W(1  -  J2)1^  q^Cl  -  W2)1^  '2(1  +  12)  (2.3) 

i  log  Tt  (1  -  l2)  <* 


(2.4) 


Although  adiabatic,  conpressible  radial  flows  are  restricted  to  be 
outside  a  Uniting  circle  (or  sphere)  on  which  1^  =  (7"  l)/(7  +  l), 

(M  =  l),  diabetic  flows  are  not  so  restricted.  For  integration  of  equa¬ 
tion  (2.3)  gives 

rn  W(1  -  B2)7'  _1  =  (  r^t'l  -  W2)1^  ~3  (1  4  **)  d  r  -»  c| 

Both  Cj  and  the  B.H.S  of  equation  (2.5)  are  not  negative.  Since 
W(1  -  w2)  xh  ~x  cannot  exceed  its  maximum  value  s 

vwVf11 


we  obtain  the  inequalities 

rr*(i  - -ar2)1^  1(i+|4r1!|2)^r+ci^c2rI 


(2.6) 


If  q^  =  0,  >“  0,  (adiabatic  case),  the  flew  is  characterised  by  an 

inner  limiting  circle  at  rn  For  arbitrary  functions  q^r)  "w® 

would  expect  one  or  more  circles  limiting  the  flow  on  either  the  out¬ 
side  or  the  inside,  (cf.  discussion  of  D  =  flows  later  in  this 

section.)  For  many  functions  q_.(r)  wo  should  find  no  limiting  circles. 
Tims  a  bounded  function  q^(r)  can  be  found  giving  HO  r  (X^  2  l)  or 
giving  (1  -  HJ2)  0(  r  2  (K^.  l)  at  the  origin.  If  W  is  to  be  differ¬ 
ent  from  zero  at  r  =0  then  there,  and  also  /  r  -*-0  if 

0  <■  •<£  1  or  if  0  <  <.  7  -  1.  Ins  order  that  no  limiting  circle  ex¬ 

clude  the  origin  it  is  i^cessary  if  q^-^Cj  r*3  near  the  origin  either 
that  K-  exceed  -1  or  that,  for  K„  —  -  1,  C  ?>7  C<i/a  [7  -  l). 
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Tc  illustrate  seas  of  these  remarks  we  shall  consider  two  types  of 
radial  flcar.  In  equations  (2.2  -  2.4)  specification  cf  nay  one  of  the 
four  quantities,  p^.  IT,  a^.,  as  a  function  of  r  should  in  principle 

permit  calculation  of  the  ethers.  Seas  combination  of  the  variables 
night  equally  -well  be  specified,  for  example  q^ 

However,  calculations  aro  especially  simple  if  U(r)  is  specified  and  q^, 
p^,  Yt  ccnputed  subsequently,  and  this  procedure  will  be  adopted  here. 

For  our  first  example  we  shall  therefore  assume  a  unipl*nar  flow  (n  =  l) 
in  which  the  fora  of  W(r)  is 


2  -1/2  2  1/2 
*  =  r(l  *  t2)  ;  (M  =  (-^r) 


(2.7) 


corresponding  to  a  flow  extending  from  r  =  0  to  r  =  oo  with  a  continuous 

2 

transition  through  the  Telocity  of  sound  at  r  -  (/  -  l)/2.  The 
heating  parameter  q^.  has,  from  equation  (2.3),  the  value  2  at  the  origin 

and  depends  anr  according  to  the  equation 

*  ft  '*]  “  [■  *  '*]  "‘A 


(bring  to  the  divergence  of  the  stream  tubes,  q^  passes  through  zero  for 
a  velocity  greater  than  that  of  sound,  namely  for  Y2  =  (y  -  l)/(7  +  l), 
(M2  =  4/(3  -  y  ),  r2  =  2(  y  l)/(3  -  y  ).  At  sonic  velocity 


q|f=[2{r+l)]’':l/2 


(2.9) 


in  accordance,  with  equation  (3-3)  of  D-l.  The  r-dopendonce  of  Vt  is 
found  by  integration  to  be 

(YtAt  jnRx)  =  (^)[z7 (  r  -  lj]1/2  r2(l  +  r2)  ‘  (7-1  +  2rr2) 

(2.10) 


-3/2 


in  which  reaches  its  maximum  valuo  whore  =  0 .(*)  The  velocity 

ratio  V/V.fc  is  obtained  by  multiplication  with  Tt.  The  other  quantities 
can  be  obtained  similarly. 

In  figure  1,  a,  b,  c  the  computod  characteristics  of  this  diabatic 
source  flow  are  illustratod  for  y  -  1.400.  The  net  energy  change  of  this 
fluid  in  passing  from  r  =  0  to  r.  =  QQ  is  moasurod  by  pT.  (<*»  )/T.  *1  = 

Vt(oo)A{.  max  =  0e708.  Thus  only  part  of  the  energy  added  in 

accelerating  the  fluid  from  rest  and  zero  temperature  at  r  =  0  to  maxi¬ 
mum  energy  at  r  =  0.707  is  removed  in  continuing  tba  oxpansion  to  where 

V  =  V  ,  U  =Oo,  R  =00.  The  maximum  rate  of  heat  evolution  occurs  for 
t 

(*) - 

Actually  V^.  could  be  uifferent  for  each  streamline,  in  analogy  with 
.the  generalized  adiabatic  flows  considered  in  references  4,5,  and  6, 
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FIG.-,  la.  RADIAL  DIABATIC  FLOW 


r 

FIG.  lb.  RADIAL  DIABATIC  FLOW 
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i 

r  —  0.36  (M  ~  0.80)  whore  it  oauois  (for  V.  =  10  unit  longths/soc) 

^2  x  ilr  see"*  ,  expressed  as  fraction  released  per  second  of  total 
heat  added  between  r  =  6  and  -r  ?  0.707.  The  value  of  IS  at  r  =  0.707 
where  Q  =  0,  is  1.58.  For  lar^or  values  of  r,  (J40  but  approaches 
zero  rapidly  as  r-wco  .  The  stagnation  pressure^  Figure  lb,  which 
reaches  a  isininarn  at  r  ■==  0.707,  increases  indefinitely  a3  r-VroP  , 
illustrating  the  extreme  sensitivity  of  hypersonic  flow  to  heating. 

Thus  our  example  shows^tho  possibility,  already  discussed  in  the  hy¬ 
draulic  'approximation,  of  obtaining  supersonic  flow  with  combustion 
even  though  the  nozzle  has  no  converging  part.  One  could  have,  . or  ox- 
oaple,  increased  I!  beyond  JL.58  by  adiabatic  expansion  instead  of  by  the 
exothermic  expansion  calculated.  Such  a  radial  -flow  could  possibiy  be 
set  up  experimentally  in  a  nozzle  consisting  of  two  parallel  discs  be- 
tween  which  combustible  mixture  is  fed  through  “tubes  on  the  axis  of  the 
discs. 


The  second  case  of  radial  flow  to  bo  considered  originates  in.  the 
use  of  potential  functions  in  two  and  three  dimensional  fields  (cf. 
Sections  1  and  3).  It  is 'instructive  there  (see  D-2)  to  study  flow*  of 
constant  type  (D  =  D  =  constant)  leading  to  a  family  of  functions 

g^(N,  Do)  parameterized  by  Dq.  Specifications  of  the  same  gjj  functions 

for  purely  radial  flow  amounts  in  equations  (2.2  -  2.4)  to  specification 
of  a  combination  of  the  variables  and  W  and  will  illustrate  again  the 

differences  between  adiabatic  and  diabatic  flow.  We  note  that  the 
second  arbitrary  function  F(^jj),  which  occurs  in  equations  (1.16,  18), 

sust  also  be  given,  if  g  is  to  be  determined.  For  the  present  simple 
illustrative  case  we  place  ?(^  ^ )  =  0  which  corresponds  (see  D-2, 

Section  3)  in  more  general  cases  to  irrotational  IT  flow  in  which  (V/il) 
can  vary  between  streamlines,  (in  the  next  section  irrotational  flows 
for  v/hich  F(^jj)  £  0  will  bo  computed.) 

We  can  again,  consider  two  and  three-dimensional  radial  diabatic 
flows  simultaneously.  From  the  earlier  theory  (D-2),  <#>  ^  must  satisfy 

the  equation 

r~n  ^  (rn  ^  )  =  0  (2.11) 


so  that 

-n/&„ 

(KA)  =  r- 

-where  K  is  the  the  integration  constant  and  N  'g  becomes  (D-2) 
N2g  -  (1  +  Do)(rr^Do  H  1^°  -  1)  D0  4  "I 

=  £•-  n  log  rj  *’1  D0  =  “1 


(2.12) 


(2.13) 


* 
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Now  for  real  flows,  H  g  must  bo  - 
r  is  correspondingly  restricted. 

,^>i 


From  equation  (3.1)  wo  see  that 


(2.14) 


for  all  values  of  Dp,  Since  the  value  of  Dp  determines  the  typo  of  flow 

in  non-radial  cases  according  to  its  sign,  we  can  carry  over  the  termin¬ 
ology  and  say  that  for 


elliptic  flows,  DQ  ^  0;  r  7>  1 
hyperbolic  flows,  Dn  <,  0;  r  •<  1 


(2.15) 


(Dq  =  0  is  not  considored  since  it  implies,  for  F  =  0,  that  streamlines 

are  parallel.)  '  Thus  a  second  departure  from  adiabatic  flow  behavior  pre¬ 
sents  itself,  for  now  we  have  exhibited  flows  which  can  be  wholly  within 
as  well  as  wholly  without  the  circle  r  =  1.  Also  where  r  =  1,  (N.^C)  =  1 

2  2 

and  N  g  =  /  l f  becomes  infinite  for  all  Dp  0)  whereas  in  the  adiabatic 
case  M  =  1  on  the  circle  limiting  the  flow  internally. 

Vortex  Flows 


*  lI=W(r,8)/S 


(E.16) 


the  equations  now  are  (cf.  (D-l),  equations  (2.5),  (3.2)  with  \/ •  =  0), 

in  cylindrical  polar  coordinates, 

—1  ^  1®G  Pa.  P  y 

r- - =  <2-17> 

b  log  p,  -  — ,  oj  vr 

■  <#5>  rrrsn  <■•”> 


3  -  V-1'  [x  -  n  J 

-*1  (1  +  $  -  0  W2)’ 1(1  -  W2), 

^  =  (y+D/(r- 1) 

b  l0S  T  z 

(2  r)  -je-—  =  (1  -  W2)  q/T 


(2.18) 


(2.19) 


(2.20) 


These  equations  sh6yi  that  for  an  adiabatic  vortex  flow,  tho  velocity, 
stagnation  pressure  and  stagnation  temperature  are  independent  of  0. 

For  an  irrotafcional  adiabatic  Y£  field,  p^  is  everywhere  tho  same,  and  Y/r 

=  constant,  thus  leading  to  an  inn^r  limiting  circle.  For  a  rotational 
adiabatic  flow  Vf  -  W(r)  ec  and  equation  (2.18)  can  be  integrated, 
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(2.21) 


The  functional  dependence  of  7  audio  ^  upon  r  is  now  not  rostricled  so 

that  flows  In  which  there  are  no  Uniting  circles  are  easy  to  construct. 

2  2 

Tar  example  with  (compare  equation  (2.7)  et.  seq.)  w  =  r  /(l  +  r  ) 
again 


and 


2v  .  2^ 


U)  jf  =  (2  +  r  )/(l  +  r  ) 


log  (pt/pt  )  j^log  (1  +  r2) 


(2.22) 

(2.23) 


Thus  a  vcrticity  that  is  everywhere  finite  can  permit  adiabatic  vortex 
flows  without  limit  circles.. 


In  diabetic  vertex  flows  of  the  typo  described  by  equation  (2.16) 
7  and  p^,  can  depend  upon  8.  However,  v/e  notice  that  the  difference  be¬ 
tween  log  p^.  and  a  certain  function  of  7 


-yzz  f  T(!  “4  ^)(1  +4  w2)  X(1  “V)  ld  W 

s-^n-Al  +4T2]  [w2  -  l]  ^ 


(2.24) 


must  be  independent  of  9  or 


log  pfc  -  rx(Vf)  +  f2(r)  (2.25). 

7e  can  summarize  the  properties  of  the  radial  and  vortex  flows  thus  r 

(i)  rotational  or  irroiatienal  adiabatio  radial  flows  havo  limit 
circlos  owing  to  the  form  of  tho  continuity  equation; 

(ii)  adiabatic  or  diabatic  irrota1 ional  7  vortex  flov/s  have  limit 
circlos  owing  to  £  he  form  of  the  irr  o  tat  I  oria.l  i  ty  condition; 

a 

(ill)  flows  without  limit  circles  can  bo  constructed  by  adding 
heat  to  radial  flows  or  vortex  flows  or  by  changing  the  vorticity 
pattern  of  vortex  flows* 


30  IRROTATIOML  N  F1C7S  IN  770-DIMB1SI0HS 

If  tho  vector  N  -  V  £g(W)  R  t!  is  irrotational  thon  its  po¬ 

tential  ^  must  satisfy  equation  (1*19)  and  the  hocking  function  io 

given  by  equation  (l015)„  Since  thoso  two  equations  involve  two  arbitrary 
functions,  g(N)  and  F \<p  tho  nature  of  N  flow  patterns,  even  though 
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irrotational,  and  the  corresponding  distribution  of  hoat  sources  can 
still  be  quite  raried.  In  order  to  illustrate  the  nature  of  thp  flews 
wa  shall  choose  cortain  functions  g(N),  F(  <$>  and  give  approximate 

solutions  of  equation  (1.19)  with  appropriate  boundary  or  initial  con¬ 
ditions. 


Elliptic  Flows 


Equation  (1.19)  is  simplified  if  D  =  1  which  implies  that  the  flow 
is  always  elliptic  (cf.  D-2)  and  that 


g  =  2(k2  -  N2) 


-1 


(3.1) 


The  equations  become,  with  the  added  assumption  F(  <p  ^)  =  <p  ^ 

q3(l  -  *Vl2)  =  i  k2  £l  -  (Y*Z)  K2/r  <f>  „  +  »S  V  •  I/*  ^  (»•*) 
V2  (k2  +  K2) *p  3  (3.3) 

4  ' 

K  uniplanar  elliptic  flow  described  by  these  equations  is  the  first  case 
of  irrotational  N  flow  we  shall  treat. 

The  constant  \i//2  can  be  removed  by  a  change  of  scales  in  the  (x,y) 
plane* 


X 1  =  *X./J*~ 

which  also  transforms  N  j= 

N  =  N,  k//T~ 


(3.4) 


v» 

(3.5) 


The  equations  for  and  <p^  become 

2(1  -  |  »12)  qK/k2  -  fl  -  (Y*Z)  N^/2  rJytjj  +  Kj3  Vx  •  S/2  Y  ( 3.6) 


and 


+  (3.V-) 

the  subscript  1  referring  to  the  new  coordinate  vector  r^„ 


The  parameter  k  can  be  removed  by  a  scale  change  in  all  the  irrotation¬ 
al  flaws  parameter ised  by  D^. 
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In  crdor  to  '•cnplefcc  the  specification  of  our  problem  ~o  adjoin  to  (3*7) 
the  boundary  eon.iitj  cn 


*P\  '*■  cci;  °) 


^  =r,  >  1 


(3.8) 


b<p\i 

*T1 


rl=1 


(3.9) 


corresponding  to  a  diabatic  flew  around  the  circle  =  1  with  pro¬ 
scribed  values  for  *£>  ^  tipon  tho  larger  circle  r^  =  r^*.  Other  boundary 

conditions  and  choices  for  the  functions  F(  „)  have  beon  discussed 

(7)  ” 

briefly  elsewhere  .  ’.To  believe  that  usofol  results  may  be  obtained  if 
b$/  is  takei.  to  vanish  on  r^  1  corresponding  to  &  flow  that  is 

unlforn  in  direction  at  'Argo  distances  from  the  cylindrical  obstacle 

r,  ^  1.  Then  vrh-. **  the  extension  to  values  of  D  other  than  1  is  made  by 
1  o 

scaling  rules  the  flew  would  romain  uniform  in  direction  at  large  dis¬ 
tances.  Also  !>r.  Hinsdale  of  these  Laboratories  has  eheem  that  F($£  ^) 
in  the  more  general  elliptic  flcnr  equation 


(3.10) 


can  bo  so  choson  that  exact  linearisation  ia  possible*  Thus  ho  makes 
tho  transformation 


dt  exp  (-  f  F (cr)  d<?~) 


(3.11) 


which  loads  to 


V2  W  =  $  0J0  »F(^>h)  exp  F(s)  ds)  (3.12) 

Tho  conditiou  for  linearity  of  (p  (ip)  is  that 


Fs(<^>^)  -  "K  =  conotant 


(3.13) 


-According  to  tho  valuo  of  A  there  are  then  three  possibilities 


U)  a  =  o5  p;^>H)  ==  -&9  +  \jf“oc  +<f>0y hfi,  (3.H) 

<])  (t|0  =  - 


D.  L,  Hicks,  nSiittple  Diabatic  Flows  of  Elliptic  Type0*  Phya.  Rev.  72, 
179(A),  (July  15,  12-17)  „ 
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(3*15) 


(ii)  A  =  &  2;  =  Octaa^f  (fa  + 

f=/?sin[o (p  Cf)  =<*-2<r-r0> 

(iil)  A  =  -O'  *-  t(<f  j,)  =  -  or  coth  jor  +?0j]l  (3.16) 

is  coshp  (^,  *i>  j]*y0:<p  op)  =  -0'  2(  T|f-  ¥o> 

# 

Far  these  three  cases  calculation  of  diabetic  floe  fields  can  be  based, 
in  principle,  wholly  upon  analytical  computations.  Further  study  of 
these  cases  night  veil  be  rewarding  since,  for  example,  exact  calculation 
of  transonic  diabetic  fields  would  be  possible.  Our  further  treatment 
here  will  be  based  upon  approximate  linearisation  of  equation  (3.7).  The 
parameter  ^  fixes  the  scale  of  If  <f>  q  is  nsall  enough,  nax 

will  be  snail  compared  to  one  and  equation  (3.7)  becomes 


V.  '  4>* 

* 

whose  solution,  for  the  boundary  condi ticns  equations 


(3.17) 

(3.8,  3.9)  is 


£>H(*yO)  * -  [kx (l)I0(r1>  +  ri(1)K0(ri)][K1(1)Io(rit)  + 

t^^oos  ©jjl0(l)  -  ^(lij^^)  -t  [*(1)  -t  KX(1)J  Vr,)]} 

|  I0(l)  -  Ijd)]  K^*)  +  |\(1)  +  Kx( l)j  I1(r1*)j|“1  (3.18) 

Calculations  have  been  made  for  r^*  =2  and  <f>  q  —  0.2864  which  insures 
that  max  =  0.203  (corresponding  to  max  1!  ~  0.242).  The  accurate 
evaluation  of  \7  °  S  in  equation  (3.6)  requires  here,  as  elsewhere,  eon- 

putation,  of  the  expression 


JlV  •  S  =liVV  -  <p  *  <P„-  2r-2^r  <fi 


r0 


(’  «*4 
-r 


^s2  f 


00 


+  r“Vr  *.f 


(3.19) 


(or  its  equivalent) *  y 

It  would  be  difficult  to  obtain  reasonable  accuracy  if  the  calculation 
of  V  0  S  were  based  upon  measurement  of  tho  divergence  of  streamlines 
drawn  by  graphical  interpolation  or  if  it  wero  based  upon  values  of  96 

N 

ftod  its  derivatives  determined  from  say  elaxation  calculations  that  might 
be  required  in  a  non-linear  problem. 
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FIG.  2b.  ELLIPTIC  DIABATIC  FLOW 


FIG.  2a.  ELLIPTIC  DIABATIC  FLOW 


In  figures  2^*,  b„  the  atreexline  pattern  and  lines  of  constant  M, 

&.  P/P0  2%  A  *ro  show,  being  approxiwately  the  'value  of  p  at 

rl*  =  2#<^  ®  to  the  cymnetry  of  the  floe  about  the  y  axis. 

It  Is  only  necessary  to  jive  these  cur  res  for  either  half-plane. 

Hyperbolic  and  Parabolic  Floss 

K**aples  of  wholly  hyperbolic  or  parabolic  floss  can  be  easily  con¬ 
structed  shore  a  Glauert  -  Prandtl  type  of  linearised  treatment  la  per¬ 
missible.  In  this  my  me  shall  obtain  sons  characteristics  of  diabatie 
fl<ars  in  curved  ducts  which  will  serve  to  introduce  an  approximate  general 
similarity  traatcent  of  ir#>taticoal  diabetic  flows •  Let  us  consider 
first  a  wholly  hyperbolic  flat  mentioned  in  D-2  {Section  8)  which  has  the 
/  proparty  that  the  flow  -velocity  is  everywhere  sonic  (M  c  l)  sad  that 

*jS  (3.20) 


This  flow  is  derived  by  choosing  g  =  for  which  D  =  -{^4  l). 

The  linesrixed  differential  equation  tor  <f>  ^  (in  two  dimensions)  be- 

eowea,  ifFty^)  (cf.  D-2,  equation  (8.4)), 


2  .  2 


An  elementary  solution  of  this  hyperbolio  equation  is 

=»0  x  coah  W-y-  »0  f) 

satisfying  the  boundary  conditions 

^,=*0 

on  y  =  0 


(3.21) 


(3.22) 


(3.23) 


Tho' linearisation  is  valid  -ten  ix  )  ~H  ~  HQ  or 

roughly  In  tho  region  )[(/■+  l)/y]  V?  »^y  j«  |!  [(r+  D/yj  1/2K02xy(^-l. 

The  expression  in  Cartosinn  coordinates  that  is  analogous  to  equation 
(3.19)  is  n 


II3  V  •  S  -  /  y2  <f>  ^  -  20x  jlSy  + 

which  in  our  present  .pproxiraation  becomes 

V-  s  =  (  y  +  i)y'Vz  no  .  2  xanh2 


^rVyy  <’'24> 
* 

Yj)  (3.25) 
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FIG*  3*  HYPERBOLIC  DIABATIC  FLOW 


in  which  the  scale  change# 


*1  " 


f  V * 
'ZtTN  y 

y  *0  y 

'  4 


is  as) 


have  been  made. 


She  heating  factor  q?J.  is  given  by 

%  =  |  K02-r  +  l)7”i/2  ^  ccah  yx  (X  -  2  tanh2  yj 
The  pressure  variation  is  found  to  be  (cf.  D-2  equation  (8* 61} 
log  (p  N^)  4  *  •-  constant 


(3*27) 


(3*23) 


fliers  fore 


*®S  (?/Pc) 


1  rv 

-  -Jv'C'Xeg  o c« h  ~  Xj  cnsJf  y^v 


?o> 


Tercperaturev  (or  velocities  &  7*/  )  can  also  bo  computed,,  The  results 
are  illustrated  in  figure  Z»  The  regions  of  validity  cf  this  solution 
are  indicated -by  the  accuracy  limit  curves  xn  ■’•he  figure  vhich.  arc  locii 
of  constant  values  of  tanh  y.,  •  Thus  at  all  oojr.ts  bonoath  tho 

0*2  curves*  the  y  component  of  JJ  is  not  greater  than  '—30$  of  the 
3t  component. 

A  more  general  analysis  of  the  hyorbolio  equation  (3*30)  can  be 
based  or.  the  Rienann  function  for  the  equation  which-  leads  to  tho  ex¬ 


pression 


Vyx 


I,  (Z ) 


<fiirSh°y v  ~  i j_f  -;':i  ■*  :fx.s  '*  r'xi  ~  '*]  +  -f  /  1-*  “"2 


cryi 


xryi 


ed  j)  i0(z)  d5  x 


(3*30) 


•where 


=  [*i- 


(•2  > 


cor.reaponding  to  the  boundary  conditions 


^  d>  v 

Ay,  Jl; 


on  y,  ~  0 


^-7  VO\ 
.  Vi  •  ll«-.  / 
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Tlu>  functions  i0(4-  ftc'  «*a  ®o3ifi-*d  rosso-*  functions.  (If  F(^g) 

=  the  solution  jafry  bo  obtained  by  replacing  I.  {2)  by  J,  .'z)J  end 

I0(«)  by  <^(2)  in  equation  (3,30),)  Tho  expression  for  $>  ,j  can  bo  trens- 
formed  to  road 


^>s(x1.y1)  =  \  [r(x!  +  ya)  +  f&q  -  yx)] 


In  this  fora  it  is  clear  how  the  nature  of  the  solution  differs  free 
.•  that  for  a  simple  wave,  Integration  of  the  values  of  =  f  as  woll  as 

of  (  <^j|/  3  jj)  =  g  along  the  support  curve  =  0  are  included  and 

the  values  are  weighted  with  the  greatest  weighting  near  tho  endpoints, 

(x^  t  y^),  of  tho  region  of  influence  for  the  point  (xj,  y^).  The 

solution  cannot  satisfy  the  linearization  condition  for  all  x  and  y  bo- 
cause  of  the  divergence  of  IQ(t)  and  Ij(s)  as  oxp  2  for  largo  z •  Tho 

heating  factor  and  pressure  satisfy  tho  same  equations  (3.20),  (3,28) 
as  before,  To  also  notice  that  change  of  K  can  be  absorbod  as  a 

scaling  factor  in  oach  of  theso  cases  of  hyperbolic  flow.  It  would  bo 
interesting  to  work  tlirough  some  case  such  as  (  d&y/  d  x)oc  (l  -  tanh  x) 

because,  since  q>T  then  need  not  be  negative  for  x  <  0,  a  greater  part  of 

the  field  would  correspond  to  flow  with  (exothermic)  combustion.  Asym¬ 
metrical  flats  would  occur  if  (  Ipy/bY)  Z"0  on  y  =  0„ 


As  a  second  example  of  duct  flew  consider  tho  wholly  parabolic  flows 
which  occurs  only  when  P  s  0,  Unlike  Uw  special  hyperbolic  flow  just 
considered,  in  this  parabolic  flow  the  Mach  nuabor  varies,  in  tho 
Glauor^-Prsxidtl  approximation,  the  partial  differential  equation  for  „ 
it.  in  thi3  case  (of,  D-2  equation  (7-3)  with  =  <f>  m)° 


(3,3*1) 


whoso  solution  is 


<fi  ^  =  f(x)  cosh  (/k  No  y)  =  k”1  V*  £ 


(3,35) 
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for  the  boindary  cpnditioos 


(3.36) 


and  is  a  valid  approximation  in  the  regions  where  j  f(i) 
jr(x)  rl  <<  {t'/k.j1/1.  The  d^ecial  solution  aquation  (3.3l)  of  the  pre¬ 
vious  hyperbolic  equation  was  of  the  same  fora.  However,  the  coordinate 
‘scale  changes  are  no m  different  and  the  sealing  of  the  heating  factor 
and  the  pressure  are  also  changed.  Thus  corresponding  to  equation  (3.35) 
and  with  the  y  scale  change  y.  =  (kl  )V^  y  wo  find  (cf*  equation  (7-5), 
D-2)  L  0 


q„  =  UjOy  tM  *  q  20ro)  f(*)J  cosh  ya  (3.J7) 

»iw«  ^(h,)  and  ^  2^*0^  aro  functions  of  (if  *>o  in 

equation  (3a37)  were  replaced  by  the  value  of  N  calculated  Aron  equation 
(3.35),  j  .nd  p/p o  might  be  calculated  with  slightly  greater  accuracy 

in  same  instances*)  The  pressure  variation  is  given  by  (cf*  D-2,  equation 
(7-6)) 

(p/po)  =  (k  -  No)  exp(-  j  <f>  VZ)  (3.38) 

For  this  parabolic  flew  it  would  be  possible  to  specify  a  function  f (x) 

for  which  q„  will  generally  be  positive*  For  a  specific  illustration, 

*  *• 

however,  wo  will  tike  f  (x)  =  x»  Then 

%  =  k1/2  N0S/2(1  -  )  xx  cosh  7l  (3.39) 

whore  x^  =  (k  Hq)^2  x.  The  streamlines,  region  of  validity  of  this 

solution,  and  values  of  Ms  (p/po),  (T, /T.  )  and  Y/V  are  shown  in 
figures  4  a,  b* 

Similarity  Rules 

A  general  linearized  treatment  of  irrotational  diabatic  flows  can 
be  givon  which  leads  to  similarity  rules  for  these  flows*  Like  moot 
linearization  treatments,  the  method  cannot  make  adequate  allowance  for 
large  perturbations  of  density  and  temperature  and  therefore  is  limited 
in  its  range  of  usefulness.  Our  development  starts  from  the  partial 
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FIG.  4b.  PARABOLIC  CM  A  BAT  1C  FLOW 
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FIG.  4a.  PARABOLIC  DIABATIC  FLOW 
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differential  equation 


^  x2  ^  y 


(3*40) 


vhioh  follows  fraa  equation  (1*18)  where  \t<fi  y/j  y|<<  j P* 
function  F(^ j)  depends  only  and  not  on  V.  Accordingly  all  effaots 
of  ohoioe  of  g(N)  and  of  free  stream  value  of  V  -  N0  oan  be  absorbed  by 


the  scale  changes. 


*lmZ  ['.*  *  iDT)]  “o'1 **  <»„  A  o) 

(3.41) 

„  2  _  f  w  2  .  1  1  2 

'  71  'L  «  iTSjJT 

(3.42) 

giving  the  equation 

2  2 

±  4-nr  *  4^?  -'«V 

d  x,  0  Yl 

(3.43) 

in  idiioh  the  +  sign  corresponds  to  the  elliptic  and  hyperbolic  cases, 
Dq  3*  or  Do<_  0.  When  -  6  (the  parabolic  case)  the  soale  change 

(3.42)  gives 

2 

1  A*  % 

d  y/  VK’ 


(3.44) 


Boundary  conditions  are  usually  representable  in  the  form 

**  +/8  (x,y)  =  A  (x,y)  (x,y) 

on  curves  J”\.  With  the  scale  changes  this  becomes 

oli^xi,:ri,Na^N  =  A 

t */*  ltwV 

where  non' 

o((x,y)  =  «1(x1py1,No);  ^  (xpy)  =/^Vxlpyl®No) 


(3*45) 


(3.46) 


(3.47) 


•1/2 


A  (x.yVA^.V  -  [d  + 


31 


Thus  whenerer  one  can  vrite  do>vn  a  solution  ol’  the  system  equations  (3*46) 
and  (3.43)  or  equation  (3.44)  for  general  values  of  the  functions 
Oi  ^  then  solutions  of  tho  system  equations  (3.40),  (3*45)  can 

be  constructed  for  any  value  oi'  and  for  any  form  of  the  function  g(No), 

•p  long  as  the  same  sign  of  D0  is  kept  (or  the  value  zero  in  the  parabolio 

case).  It  is  most  likely  that  goneral  solutions  of  the  system 
(3,46)  and  (3,43)  or  (3,44)  mill  be  available  if  the  system  is  completely 
linearized.  Suppose  for  example  that  =  0^  previous  dis¬ 

cussion  of  various  cases  shows  that  three  basic  solutions  for  example, 
equation  (3,18)  in  the  elliptic  case,  equation  (3.30)  in  the  i^perbolio 
oase  and  the  expression 


85  t(x)  cosh  y1  +  g(x)  sinh  yI 

(  ^  »  f  U),  -  g(x),  <«  y  =  °) 


(3.48) 


in  the  parabolic  case  permit  iroediate  calculation  of  the  flow  pattern, 
Pto,  for  all  values  of  N^^just  by  making  the  scale  changes,  equations 

(3.41)  and  (3042)a  Thus  all  elliptic  and  hyperbolic  flcwro  (withF(^^) 

—  <^>  jj)  and  with  the  boundary  conditions  given  in  equations  (3.8),  (3,9) 

and  (3.32)  respectively  cen  be  derived,  no  matter  what  the  function 
g(No),  from  equations  (3oX8)  and  (3.30)  in  the  Glauert-Prandtl  approxi¬ 
mation.  The  equation  for  q,?  i 0 

%  ~«1  4  %  ["o-eO'o)]  <3-49> 

in  which  q,,  qr,  are  functions  of  N  that  are  detorminable  once  g(K  )  is 
•  i  o  o 

given, 

TTo  thus  find  the  following  similarity  laws  for  linearized  irro- 
tational  diabatic  flow  (in  addition  to  equations  (3.41,  3,42)) 


B(N0*)/«(K0)  =  f.\  '  s(Ho*)/H0 

w.  A 

2  1  V2 

(3,50) 

r  c?<4>n  /  jdti 

w 

(3,51) 

d<Pw  1  . r  o'?  1;  1  W 

WIJ  /lJy'1WJ\  := 

tm  *• 

(3.52) 

e(V»  ^kVpO'o-  **»>  •- 

(3.53) 

>2. 


where  (cf.  D-2  equation  (5-8),  H,(N)  =  exp 

«•# 


/ e( *)  i  *ZA. 


and 


viv-i  b2t^Jni!± 


(3.54) 


Equation  (3*53)  could  be  used  to  calculate  pressure  farces  on  bodies  in 
diabetic  flow*  The  transformation  of  Nq  and  g  would  of  course  distort 

the  body  being  considered  just  as  in  Glauert-Prandtl  transformations  for 
adiabatio  flow*  He  note  tljat  for  all  irrotation&l  diabetic  floes 

(1  ♦  %y-  B2  g)  »  •  V  log  *t  =  2%  <S-S6) 

If  is  the  free-streaa  value  of  the  stagnation  temperature  then 

log  (*tAto)  =  fa  d  ^>h/42(1  +  E2  g)  3  (3.56) 

where  the  integral  is  takan  along  a  streamline.  The  similarity  law  here, 
in  the  Glauert-Prandtl  approximation  is 

[  Jx  **  W>/^  l0S  Vo>}  =  o4(H0'...)Ho(1  4  g(»0))/ 

%W(1  +  TT  V2  8<V»  <3-‘7> 

Once  N,  p,  T^  are  known  then  of  course  M,  V,  p^,  T,  and  Q  can  be  com¬ 
puted.  He  furthor  note  that  the  vorticity  in  V  language]  V  x  J,  I 
two  dimensional  irrotational  flow,  simply  equal  to  (-  —  d  log  g  T for)'  )• 

The  aocuraoy  of  the  various  approximate  expressions  we  have  given  may  be 
expected  to  be  different  for  the  different  flow  patterns  and  algo  for 
the  various  equations,, 


4,.  PEREfitBATJON  OP  UNIFORM  FLOW  BY  A  LOCAL  SOUBCE.  OP  UEAT^ 

Among  the.  siagpior  rotational  diahatic  flows  only  one  has  boen 
troatod  fully  by  us,  that  of  an  almost  uniform  flow  deflected  by  one 

B.  Lo  Hicks,  "Perturbation  of  Steady  Uniform  Flow  by  Localized 
Sources  of  Heat",  Phya„  <ev,,  73,  636  (L)  (March  15,  1348). 
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w  CX  2 

local  heat  floored  of  the  far*  Q  C#  rl  „  More  gonoral  rotational 
problems  can  bo  discussed  by  examining  per  tur bation  duo  to  hoat  source* 
In  non-uni  for*  flows,  by  superposition  of  more  titan  one  heat  source  in 
almost  uniform  flow,  or  by  relaxation  or  characteristic  paramotor 
methods  for  high  Mach  number*  (cf .  basic  equations  -  Section  l)a  We 
hope  to  report  some  of  these  advanced  calculations  at  a  later  time* 

W»  had  expressed  a  doubt  in  the  Uadison  paper  that  a  unique  ele¬ 
mentary  source  of  heat,  analogous  to  a  fluid  source,  oould  be  defined 
for  diabatio  flow*  Comparison  of  the.  results  to  be  discussed' in  this 
Section  with  the  calculations  of  reference  (9)  lead  to  the  conclusion 
that  difficulties  such  as  infinite,  jump  in  entlvalpy  or  appearanoe  of 
▼ertex  sheets  could  be  ignored  in  a  first  treatment  if  the  total  heat 
per  second  added  to  the  fluid  is  kept  constant  as  the  spatial  extent 
of  the  heat  source  is  reduced  to  make  it  a  line  or  point  sourco.  (its 
are  indebted  to  Professor  Tsien  for  sending  us  a  prepublication  copy 
of  his  abstract*)  It  *ust  be  reoognlsod,  however,  that  the  infinite 
velocities  connected  with  an  elementary  line  souroe  make  it  partially 
unsuitable  for  a  perturbation  theory,  He  will  therefore  present  here 
the  perturbation  calculation  for  the  smoothed  source  exp(~  dr^). 

The  unperturbed  flow  is  uniform  and  of  velocity  in  the  direct¬ 
ion  of  the  positive  x  axis.  In  the  perturbed  flow  let 


V  =V  i  +  Vs  «  (V  4  u»)  i  4  v»  j 

/w  O  r*J  '  O  ^ 


PSP  +  P°  * 


/,  =  A>+/’ 


Tt  =  Tto  +  V 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


in  which  the  subscript  zero  refers  to  the  unperturbed  state  and  the 
prise  to  the  perturbed  state**  Each  of  the  quantities  V0,  po„  p  is 

constant  while,  in  the  case  considered  heroc  V®,  p®  and  p  1  are  func¬ 
tions  of  (x„  y)«  Tho  first  order  perturbation  equations  are  derived 
from  equations  (l.l)  to  (1.3).  Thus  =* 

^  The  effects  of  a  line  source  of  heat  have  boon  described  further  in 
recent  first  order  perturbation  thoory  (9,  10)„ 

/  «  e 

iy^H.  S.  Tsien  and  U,  Beilock,  HUeat  Sourco  in  a  Uniform  FlcwM,  J.  Aero, 

So*  16,  756,  (Doc„b  1 949) o 

^^B.L.  Hicks,  wAn  Extension  of  the  Theory  of  Diabatic  Flow”,  Phya,  Rev. 

77,  286  (b)  ( Jan0  15,  I960), 

For  this  first  elemontary  but  fundamental  computation  it  sooraod  appropri 
ato  to  use  tho  Y,  language  as  boing  more  familiar.  For  more  complex 
flows  or  in  more  accurate  treatments,  the  conclusions  in  D-2  (Section  2) 
suggest  that  tho  W  language  may  bo  more  useful, 

t 
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Vi*'  ‘  /%  V, 

V»  v*  v  A 

^  O 


o 

p 


r  * 

v  *>  y 

O  ^  X 


j  VJ 

(4.&) 

*• 

•  - - •  0 

(4*6) 

a  ** 

-  a- 

(4.7) 

Using  equations  (4,6)  and  (4*6),  aquation  (4.7)  osn  bo  re-written  to 
first  ardor  toms  as 


a  -  «o2) 


<* 


=  ( r  -  1)  «/  qAo2 


It  can  bo  shown  that  lO  =  J  V  *  V*  I  *•  °*  Bd0an<1  ardor*  Wo  oun 
therefore  lntorduce  ft  potential  <p  ^  for  w^ieh  ¥•  =  Vq\7^y  and  after 
asking  the  scale  changes  x  =  x^j  1  -  UQ^  f  ,  y  =  y^  (Mft  /  l)  arrive 
at  the  basic  equations  for  this  type  of  flow 


+ 


“  (r  -  i)  u02  «A03 


(4.9) 


The  notion  of  tlio  heat  aouroe  as  an  effective  fluid  sauroe  (of*.  D-l)  is 
hare  illustrated.  We  note  that  the  R,.H»S»  of  equation  (4 ,9)  is  simply 
2<l|l*  in  tho  notation  of  D-2,  evaluated  in  the  IVoo  stream  where  the 

velocity  of  sound  is  ao<> 


Lot  ue  trout  a  subsonic  case  («  <l)  and  take  Q  to  be  a  fuuotlon 


only  of  r^,  namely 

5 

Q  »  2  a  q..  a  e 
^1  o 


-or  r. 


(4.10) 


wr: 


or  aa  uood  In  the  Wisconsin  Symposium  paper  Q  ~  ?.  &  q  V.  e  • 

°  ;  o 

where  oi  ia  a  parameter  of  compactness  and  q^  moauuiuu  the  iutoiuiity  of 
the  heat  roloueo,,* 


Wo  oon  show*  if  the  perturbation  velocity  u*  is  infill  comparod  to 
V^,  how  q.^  con  bo  calculated  and  given  a  practical  moaning.  When  u*  is 

small*,  t lie  Integral  in  equation  (4.21)  can  ho  rop'luoud  b y.J*  (Q/V0)  dx 

whioh  in  easily  evaluated,,  For  larger  u®  the  siuiio  dovol opmout  still 
givoa  a  qualitative  idea  of  tho  moaning  of  <]^u 


From  oijuutlon  (4*22)#  the  enthalpy  incroauo  from  x  -  "'00  to 
x  »  4  ©o  along  tho  utroumline  y  -  0  iu 

«v/\Tt  *  oJ(Ar  -  opT0  •  2(/-l)(IT'*) /2{l  -«„*)  Vo-1qi  (4-1) 


3‘) 


Footnote  ccnticwd  fr  oat  ;«j 

If  =  (4  T^/T^)  is  used  as  a  csuure  of  the  heating  talua  of  the  com¬ 
bustible  mixture  then  in  terms  of 

,x  =  [uyar1^  r-  ua  -  *02)lyl?k10<  e-u) 

end  is  thus  directly  proportion*!,  to  k^  end  to  the  scale  defined  by  the 

compactness  parameter  of  •  The  heat  output  of  a  combustion  chamber  is 
usually  expressed  in  terms  of  the  loading  factor  0  9  the  rate  of  energy 
release  per  unit  Toltaee  end  pressure  in  the  ch&nber.  If  tQ  is  the  time 

necessary  for  combustion  to  bo  completed,  then  it  my  be  sham  that 

8  =  [r/(r-  u]  (k,A0)  u-iii) 

or  if  /S  *  is  the  leading,  given  in  the  usual  practical  unit 
BTU /f t^-atm-hr  and  t^  is  in  seconds. 

/$  •  *  3.44  x  104  (k1/to)  (4-ir) 

The  two  measures  of  heat  release,  q^  and^  ,  can  be  related  if  we  say 


that  the  distance  A  x  =  (V  t  )  necessary  to  complete  combustion  is 

-1/2  0  0 

kgOC  7  where  is  of  the  order  of  two.  He  then  find 

«i  =h  V'2?>1/2/(1 "  “0?)1/2J  W  <4-’) 

*■*  ^/ji 

For  q  in  ft,  CK  in  ft*  , *  in  the  practical  units  and  with  — 2, 
-l/2 

V  (1  -  M  2)  =  0.130,  /=  1.40,  V  =  150  ft/sec,  we  obtain 


q^'*''  3.5  x  lO^Vd) 


(4-ri) 


Therefore,  for  a  combustion  chambor  for  a  turbo-jet  engine  for  example 

with  /$*=■  2  x  106  BTU/ft3-at»-hr,  (X  =  0.44  (  A  x  =  3  ft),  qx  is  0.16, 

which  corresponds  to  (l/6)  of  the  stoichiometric  value  of  0.96  cal- 

o  **1  ^2 

culated  from  (4~ii)  for  the  same  values  of  MQ(l  -  )  ,  Y  CX  . 

Equivalent  results  are  obtained  if  the  volume  of  the  combustion  space 
rather  than  its  longitudinal  dimension  Ax  is  connected  with  cX  . 

The  value  CX  s=  1  i#  used  in  the  plots  to  be  discussed  later*  If  a 
different  value »of  C*  is  of  interest,  tha  x  and  y  scales  must  be  ro- 

plaoed  by  CX  */2  x  an<i  &  */2y,  by  (q^  &  ),  T*  by  (y/S~T *)p/0*  bjr 

(yR/5*)  And  u®  by  uv)n 
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and  thickness  (the  Q  of 

1/2 


The  total  heat  added  to  the  fluid  per  unit  tii 
reference  (9))  is  thus 

^  jy2  00 

„s«  **  «  (1-*o2>  p  «  rXdrl  =  **(».V  V.>«5 

S  (4.11) 


=  ZT// 


ML 


0  o 


Ao  oar  ding  to  equation  (4.10),  q.  la  a  function  of  r,  «ly  and  equa¬ 
tion  (4.9)  reduces  to 

rf  1  (rx  =  2(y-l)  a*;1  qj  .‘<rri  (4.12) 


Okie  integration  gives 


^  =  (r-o.0,1r1-W‘S 


(4.13) 


which  differ*  fra.  the  usual  (fluid)  point  aowea  in  two— diaona ions  in 
that  there  is  no  singularity  at  r^  =  0.  As  a  consequence  it  would  he 

possible  to  insure  the  "validity  of  the  perturbation  treatment  in  all 
parts  of  the  flow  for  not  too  large  values  of  q^.  It  is  iwportsat  to 

note  that  S/fi  with  respect  to  the  space  has  radial  components  only, 

hut  that  the  y  <fi  with  respeot  to  the  physical  space  is  the  perturbation 
velocity  vector  which  is  not  syssaetrical  about  the  heat  source  location 
in  either  the  physical  or  the  r^  space.  However,  for  Uach  numbers  which 

are  very  close  to  sero  the  perturbation  velocity  deviation  Aron  syanetry 
in  both  the  physical  and  space  is  very  snail. 

Because  of  tlio  irrota^onality  of  V*,  integration  of  equation  (4.5) 
for  t)»  perturbation  pressure  p*  is  imediate,  and  the  vanishing  of  u* 
at  infinity  yields 


p,=*A  To  u' 


(4.14) 


This  simply  expresses  the  pressure  variation  needed  to  effect  the  change 
in  momentum  in  the  x  direction  associated  with  the  heating,  other  com¬ 
ponents  of  the  momentum  being  negligible.  Finally,  elimination  of  ^  •  VJ 
between  equations  (4.6)  and  (4.8)  followed  by  integration  loads  to 


/** 

- —  s=  U  a  u‘ 

Pa  °  0 


(4.15) 


+  (iroc)1/2(y-i)(i-Mo2)1/2  m^1  qi «  7l  [1  +  °r-(a  1/,z  ^ 
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where  P  •  — 0  as  1^-^-  -»  .  The  expressions  for  the  perturbed 
quantities  in  the  .case  of  a  line  source  are  obtained  by  letting  C*— >-oo. 

If  =  0.20,  Of  —  i.o,  y  =  1*40,  and  =  0. 098,  (corresponding 

to  kj  =  (Alt/fo)  =  Oi. 68)  then  aex  j  (Y'/V^)  |<(  0*13,  and  to  this  extent 

the  linearization  is  justified.  Stroaalines  for  this  case  (if  2  neglected 
o capered  to  unity)  are  shown  in  figure  5a  together  with  lines  of  con¬ 
stant  Telocity,  (or  of  perturbed  statio  pressure)  and  also  lines  of  con¬ 
stant  perturbed  temperature  (either  static  or  stagnation  within  the 
accuracy  of  the  perturbation  treatment). 

Owing  to  the  low  Mich  xtusber  the  static  and  stagnation  pressure  changes 
are  amall.  In  figure  5b  tho  variation  of  these  pressure  changes,  P*/p0 
»nd  p'^/py  along  the  x  axis  are  shown.  Perhaps  the  most  interesting 

aerodynamic  effect  illustrated  is  the  slowing  down  of  the  fluid  upstream 
of  the  heated  region  as  though  the  fluid  were  approaching  a  real  obstacle. 
Far  larger  values  of  q.  there  would  be  a  stagnation  point  and  actual  re¬ 
versal  of  the  flow. 


The  calculation  of  p*^  for  the  Madison  paper  was  based  on  the  simple 
formula 


£-)  -  <?b 


(4.16) 


and  is  the  curve  I  for  =  0.098  in  figure  5b.  Since  recovery  of  stag¬ 
nation  pressure  is  impossible  so  long  as  no  heat  is  abstracted  from  the 
flow,  the  apparent  recovery  must  arise  from  some  inaccuracy  in  the  cal¬ 
culation.  In  order  to  chock  this  point  a  more  careful  calculation  of 
p*^  was  made  by  integration  numerically  oquation  (1.9)  which,  for  the 

present  case,  reducec  to 


in  which  W  and  (T^/t)  can  be  computed  IVom  the  formulae  already  given 


for  tho  velocity,  pressuro  and  density  variations.  Values  of  p^.  ob¬ 
tained  in  this  way  are  plotted  as  Curve  II  (for  q^  •=  0.098)  in  figuro  5b. 

Comparison  of  1  and  II  shows  tho  large  inaccuracy  in  oquation  (4.16), 
which  is  undoubtedly  duo  to  the  large  density  variation  along  tho  x  axis, 
(of.  figure  5e0 )  As  a  furthor  chock  the  calculations  wore  repeated  for 
q^  —  OiCO.ll 5  which  is  small  enough  so  that  equation  (4.16)  should  bo 

reasonably  accurate.  The  p8^  ourves  for  this  value  of  q^  for  tho  two 
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FIG.  5b- ALMOST  UNIFORM  FLOW 
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methods  of  calculation,  ns  seen  from  figure  5b,  do  indeed  lie  olose  to¬ 
gether*  From  equation  (4*16)  the  asymptotic  expression  (x  — >-+  oo  .  is 


hr-)  =  -r  ( ?-nfiF  *A1  -  M„z)l/Z< 


(4.18) 


o  aayjnp 


A  mere  accurate  asymptotic  expression  can  be  derived*  The  definition  of 
p^  together  with  the  fact  that  p  — wpQ  •«  x— < w°°  gives 


<?-> 

asymp 


/  1 


-  J)  Y/r  X-(i  -•  2^  1  <4-19> 


Since  also  V  V  as  x  — 

o 


(4.20) 


“oc  =  V  <  W 

From  the  energy  equation  (l*3) 

[~ «A>  d*  =  VTt «  ~  Tt.}  =  VT~  ~  (4-21) 

Making  the  approximation  that  V  in  the  integral  we  find  that 


lO0 


T 


=  l  +  2(7 -x)0  /a  -  »02  qx  V1 


(4.22) 


where  the  denominator  is 


‘  -  <£* 


asymp 


obtainable  from  equation 


(4*15)*  Combination  of  equations  (4.19p  20p  22)  pormits  calculation  of 
values  of  (p^f/po)  which  are  in  error  only  because  of  the  small  difference 


between  V  and  V  0 
o 

expression 


Ylhon  1,  this  combination  leads  to  the  simple 


£') 

P 

ro  asymp 


y  u 2  (i  ~  :2_) 

*  o  u  Too7 


(4*23) 


which  reduces  to  equation  (4*18)  when  T9/Tq  and  <£.< 


The  variation  of  temperature  or  density  along  streamlines  is  approxi¬ 
mately  as  shown  in  figure  5cp  the  letters  au  b„  cp  dp  on  the  curvec  re¬ 
ferring  to  streamlines  in  figure  5a„  The  donsity  or  temperature  ohangeo 
along  the  axis  amount  to  7($  of  the  unperturbed  values  and  the  correspond¬ 
ing  entropy  incroase  is  about  1/2  c^*  Sinco  max  £  log  T^/ )  s  ~  Zq^Y^Q 

the  maximum  fractional  time  rate  of  change  of  energy  of  a  fluid,  particle 
is  (for  V.  -  500  unit  longths/sec)  equal  to  about  200  per  sec0  Thus 

much  lower  burning  rates  than  were  involved  in  the  special  caso  of  racjial 
flow  of  Section  *2  produce  appreciable  perturbations  of  the  uniform  fields 
of  flow* 


t.2 


In  order  to  Illustrate  the  effect  of  Kach  nunber  let  us  take 
0.80  instead  of  0.20.  Both  the  heat  source  function  and  the  per- 

will  be  distorted  slnoe  r ^  =? 


turbation  -velocity,  bfting  /unctions  of  r-,, 
r  9\  “1  2  21 * 

I  (l  -  II  )  x  +  y  J  •  The  magnitude  of  the  apparent  fluid  source 

is  reduced  equation  (4.12)  by  the  factor  .2  j/l-(«2)^/.8  y/l-(.8]^  -  0*41 
and'  the  maximum  Telocity  perturbation  is  increased  by  a  factor  of 
(l-.04)/(l-*64)  =  2*7  for  a  given  total  amount  of  heat  added  to  the 
fluid  per  unit  time.  The  accuracy  of  the  perturbation  calculation  as  a 


whole  is  thus  about  the  same  for  =  0*8  as  for  Mq 
•values  of  q^. 


0*2  for  nail 


5*  OW  THE  COUPUTATIOH  OK  GEtfrRAL  UNIFLATIAR  DIABATIC  FLOWS 


The  mathematical  struoture  of  the  basic  diabetic  flow  equation  In 
W  language,  equations  (l*9  -  l.ll)  has  immediate  physical  implications. 
If  we  regard  W  as  the  quantity  to  be  determined,  then  we  see  that  eq¬ 
uation  (1.10)  looks  like  an  equation  of  continuity  for  a  medium 
possessing  apparent  fluid  source  proportional  to  the  heating  faotor  q^* 
The  equation  of  motion,  equation  (1.9),  however,  shows  that  the  only 
important  dynamical  quantity  is  the  stagnation  pressure  p^_  which  acts 


as  a  potential  for  the  difference  botwoon  vorticity  and  heating  terms. 
In  fact  variations  of  p^  along  streamlines  are  associated  only  with 

heating  (between  shock  fronts)  and  variations  perpendicular  to  stream¬ 
lines  with  the  distribution  of  the  (reduced)  vorticity  L 

The  vorticity  and  heating  effects  are  actually  coupled,  both  beoause  of 
the  continuity  equation  and  because  p^.  acts  as  a  potential  for  them. 

This  can  be  seen  in  dotai£  after  is  eliminated  by  taking  the  curl  of 


equation  (1*9).  Since  we  shall  only  concern  ourselves  here  with  uni- 
planar  flow,  the  curl  of  equation  (l,,9)i  yields  the  one  equation  (cf . 


D-l) 

2(1  +^Tw2)%t>w+‘Jir(u'^r 


(X  -  B2)  q„(u -j-  -  T  -f-)  log  ^ 


(S.l) 


whore  u  and  v  are  the  components  of  the  voctor  IT.*  But  q^.  is  known 

from  the  equation  of  continuity  in  terms  of  If*  Therofore  wo  can  obtain 
finally  an  equation  in  u,  v  and  their  derivatives  only  • 


-v(i-  y~  u2-v2)uxx+2u(l+  ~~  v2)Uxy+  v(1+u24  ptl  v2)  u^ 
-u(l+  ~~~  u2-4v2)vxx-2v(1+  u2)Vxy+  u(l-u2-  v*)  Vyy 


(5.2) 

+  S=  0 


Noto  that  u  and  v  here  bear  no  direct  rolation  to  the  u‘,  v1  of 
Soction  4. 
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where  .£2-  contains  no  second  partials0  Tlilc  equation  must  be  satisfied 

by  the  components  u,  v  of  the  Crocco  teeter  Vf  for  all  stoady  uniplanar 
diabetic  flows.  Functions  uyX^y),  v(x,y)  which  satisfy  it  will  determine 
the  nature  of  the  heat  source  function  q^(xPy)  through  the  continuity 

equation*  Re  can  then  formulate  a  general  method  of  attacking  diabatic 
flew  problems  if*  as  in  the  case  of  irrotational  diabatic  flows*  we  first 
specify  partially  the  flow  pattern  W  and  then  calculate  q_^>  (ll,12,13)*e 

It  is  recalled  that  specification  of  q^  first  may  lead  to  such  undesired 

coc^lications  as  slip  lines  or  discontinuous  temperature  jumps  in  the 
flow  but  that  (cf*  Introduction)  calculation  cf  an  indirect  problem  is 
generally  easier  than  of  the  direct  problem. 

Let  us  suppose  that  u  is  given  as  a  function  of  (x*y)  throughout 
the  region  R  and  let  us  inquire  what  is  involved  in  finding  v(x?y)* 
Equation  (5-2)  is  then  a  second  ordor  partial  differential  equation 
for  ▼  that  may  be  shown  to  bo  always  hyperbolic.  Consequently  there 
will  be  no  difficulty  with  change  of  typo  of  tho  equation  either  in  the 
field  of  flow  or  in  going  from  one  problem  to  ariothor.  Such  equations 
of  constant  type  occurred  as  special  cases  in  irrotational  flew  (cf. 
Section  3).  One  needs  to  specxfy~v  and  its  normal  derivative  along 
curves  j  '  in  R  in  order  to  permit  solution  of  the  equation  for  v  with- 

in  some  6ub-region  of  R  that  is  partially  bordered  by  I  Wo  note  that 

the  constant  hyperbolic  type  of  equation  (5-2)  is  not  an  obvious  advantage 
in  the  adiabatic  case  (i.e.*  where  q(x,y)  is  specified  to  be  S.  0),  be¬ 
cause  in  genera]  it  is  there  still  necessury  to  satisfy  a  second  equation 
in  v  which  dees  change  type*  However,,  it  may  be  that  some  liXtoarizod 
adiabatic  problems  could  best  bo  handled  by  finding  a  solution  of  a  re¬ 
lated  diabatic  flow  problem  containing  parameters  which  are  v&ried  in 
such  a  way  that  max  <  1  or  f\^\  d  x  dy  or  some  other  expression  of  the 

heating  factors®  importance  is  minimized,, 

V/o  can  examine  tho  equation  (5-2)  further  in-  the  limiting  cases  of 
small  W  ("incompressible”  flow)  and  of  tho  almost  uniform  flow  approxi¬ 
mation. 

b.L.  Hicks*  "On  tho  Calculation  of  Stoady  Diabatic  Plows"*  Phys,  Rev. 
74,  1230(A)  (Nov.  1940). 

VeP.  Starr*  !}A  Mathematical  Thoory  of  Convection"*  J«  13a  to  or*  6p 
188-192  (Juno*  1949). 

G*Yf*  Platzmann*  "An  Example  of  a  General  Integral  of  the  Hydro- 
dynamical  Equations"*  Private  communication. 

*  v 

Professor  Starr  has 'shown  that  initial  specification  of  tho  momentum 
distribution  leads  to  a  partial  differential  equation  for  tho 
specific  volume  which  Is  second* order  but  aly;ays_h^porbolic  and 
linear  and  calculates  a  specific  example,.  Professor  Plutzmimn  has 
goner aii 2 od  theao  results.  We  are  indebted  to  those  authors  for 
sending  us  pro-publicfttion  copies  of  thoir  papers. 


For  \T  «.  1,  the  equation  bocumes 


(-  vu  +  2uu  +  vu  )  -t  (-  US'  -  ?.vv  +  uv  )  =  0 

xx  xy  yy  xx  xy  yy 


(5.3) 


Supposing  again  that  u(Xf,y)  is  known,  tho  characteristics  of  this  partial 
differential  equation  always  exist  sinco  it  is  hyperbolic,-.  _Tfkr  -twg 
families  of  characteristics  intersect  orthogonally.  Non-liaeprit^  enters 
only  owing  to  the  coefficient  v  of  tho  tor  m  v  *  In  our  earlier  formula- 

tion  of  tho  theory  of  lew-spood  uniplanar  diabatic  flow  (D-l,  equation 
(l.  16) )  the  nonlinearity  entered  in  tho  more  complicated  coupling  equajtiohf 

ijT1  a,  o.j&J  ; 

It  might  be  instructive  to  take  computed  values  of  u(x,y)  and  v(x,  o)  from 
the  flow  of  .Section  4  and  recalculate  v,  q  '  Q,  p,  etc.,  throughout  tho 

flow  field  as  an  improvement;  of  the  linearized  treatment.  In  this  check, 
the  coefficient  v  of  v  could  bo  initially  assumed  to  have  the  value 

calculated  in  Section  4  in  starting  an  iterative  solution  for  v(x,y). 

Less  ambitiously,  the  terras  in  equation  (5.1)  could  be  computed  and  tlioir 
sum  compared  to  zero  as  a  chock  of  the  linearized  treatment  of  Section  4. 

For  flow  which  is  almost  uniform,  u  --  +  u*,.  v  -  v’,  J  u’J  ,  |v'|«u0‘ 

5=  const.,  a  scale  change  loads  to  tho  simpler  equation 


X1X1 


Vi 


with  x 


c  /,  ,  /  -r.it  c  \  o  c> 

i  ^  U  +  y-r«6  )  x  ■"  y3  -  (1  ~  uo  > 


zC1  2 


(5.4) 


The  function 


is  a  linear  combination  of  the  form 


=  AX(uo>u'xiu’yJ  +  V^^x/'x  4  A3(uo)u’y1T,y1 

y 

>f  A .  (u  )vu  v « 

4  oJ  y-j 


(5.5) 


which  normally  would  bo  small  compared  to  individual  terms  in  tlio  L.  II.  S. 
of  equation  f5.4)n  Vo  note  that  in  tins  Glauort-Prnndbl  approximation 
tho  partial  differential  equation  for  v®  has  boon  wholly  linearized  ox- 
cept  for  the  prosence  of  the  (small)  term  in  v®  vr  . 

X1  yl 


It  appears  that  bore,  as  for  irrotational  flows,  tlie  almost  uni 
form  or  Glauort-Prandtl  type  of  treatment  will  lead  most  quickly  to 
approximate  calculations  of  interesting  flows.  Vo  note  that  in  tho 


o*s©  of  all  thro©  o qua t ions  \5,2,  3,4)  the  solution  n«>utl  not  bo  ob~ 
tainod  to  at  high  a  dogree  of  accuracy  as  must  a  potential  funotion,  for 
tho  latter  Mist  be  differentiated  in  order  to  plvo  -valuos  for  u#  T,  p# 

T,  etc,,  -which  are  the  primary  -variables  of  interest. 
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